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Abstract 

This  paper  introduces  an  implicit  framework  for  the  analysis  of  uncertain  systems,  of 
which  the  general  properties  were  described  in  Part  I.  In  Part  II,  the  theory  is  specialized 
to  problems  which  admit  a  finite  dimensional  formulation.  A  constant  matrix  version 
of  implicit  analysis  is  presented,  leading  to  a  generalization  of  the  structured  singular 
value  fi  as  the  stability  measure;  upper  bounds  are  developed  and  analyzed  in  detail. 
An  application  of  this  framework  results  in  a  practical  method  for  robust  7^2  analysis: 
computing  robust  performance  in  the  presence  of  norm-bounded  perturbations  and  white- 
noise  disturbances. 


1  Introduction 

Part  I  of  this  paper  introduced  a  framework  for  analysis  of  uncertain  systems  in  implicit  form, 
combining  the  behavioral  approach  to  system  theory  [24],  the  Linear  Fractional  Transforma¬ 
tion  (LFT)  paradigm  for  uncertainty  descriptions  [14],  and  the  Integral  Quadratic  Constraint 
(IQC  [26,  11])  formulation.  We  summarize  the  main  ideas  of  this  formulation.  The  same 
notational  conventions  apply. 

An  implicit  system  is  described  by  equations  of  the  form  Gw  =  0,  where  w  is  a  vector  of 
signal  variables,  and  G  is  an  operator  in  signal  space.  In  this  paper  we  only  consider  linear 
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systems  (i.e.,  with  G  linear)  in  discrete  time.  The  concept  of  stability  in  this  formulation 
is  characterized  by  the  left  invertibility  of  the  operator  G.  Two  versions  were  considered: 
“Is-stability”,  meaning  that  G  is  left  invertible  as  an  operator  on  I2',  “stability”,  meaning  that 
G  has  a  causal,  finite  gain  left  inverse  in  the  corresponding  extended  space  he- 

Uncertain  systems  are  characterized  by  a  parameterized  equation  operator  G(A),  where 
A  is  a  structured  uncertainty  operator  of  the  form 

A  =  diag  [SiG^  ,-■■■,  ^lItl  •>  •  ■  •  ?  (1) 

Robust  (h)  stability  means  that  the  G  is  (h)  stable  for  all  structured  A  in  the  normalized 
ball  of  uncertainty  Ba-  A  parameterization  G(A)  of  Linear  Fractional  Transformation  (LFT) 
form  was  considered;  it  was  shown  how  the  general  case  reduces  to  the  canonical  version 

'  I- AA 
C 

In  addition  to  including  the  standard  theory  as  a  special  case,  it  was  shown  that  the  implicit 
formulation  allows  for  the  representation  of  IQCs,  which  can  be  used  to  describe  properties 
of  signals,  components,  or  mathematical  restrictions  in  a  problem  under  consideration. 

In  Part  II  we  take  (2)  as  a  starting  point,  and  move  closer  to  the  computational  aspects  of 
implicit  analysis  questions  by  considering  problems  which  can  be  cast  in  a  finite  dimensional 
setting,  in  terms  of  constant  matrices.  Different  cases  which  result  in  constant  matrix  analysis 
are  reviewed  in  Section  2.  These  include  some  special  instances  of  (2),  as  well  as  a  problem 
from  an  entirely  different  origin,  related  to  model  validation  [22,  13]. 

The  constant  matrix  formulation  leads  naturally  to  an  extended  version  of  the  structured 
singular  value  /j,  [5,  14],  which  is  introduced  in  Section  3.  The  issue  of  upper  bounds  to  the 
structured  singular  value  for  implicit  systems  is  extensively  considered;  as  in  the  standard 
case,  these  bounds  are  attractive  computationally  since  they  reduce  to  a  convex  feasibility 
problem.  Conditions  under  which  this  bound  is  exact  for  the  implicit  case  are  discussed. 

In  Section  4,  the  case  of  state  space  implicit  descriptions  is  addressed,  and  the  analysis 
conditions  are  related  to  the  case  of  linear  time-varying  (LTV)  structured  perturbations, 
already  considered  in  Part  1. 


z  =  0 


(2) 
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The  paper  concludes  with  an  application  of  the  machinery  to  robust  performance  tests  in 
the  presence  of  structured  uncertainty  and  white  noise  disturbances,  which  has  been  referred 
to  as  the  Robust  H2  performance  problem  [15,  23].  The  main  idea  is  to  consider  deterministic 
descriptions  [16, 17]  of  white  signals  obtained  by  constraints  in  I2  space.  The  method  described 
in  Part  I  allows  for  these  constraints  to  be  included  in  an  implicit  analysis  problem  of  the 
type  (2).  Some  examples  are  included  for  illustration. 

Preliminary  versions  of  these  results  were  presented  in  the  conference  papers  [18,  19].  The 
proofs  are  collected  in  the  Appendix. 


2  Motivation  for  Constant  Matrix  Analysis 


In  many  important  cases,  robustness  analysis  can  be  conducted  in  a  constant  matrix  repre¬ 
sentation.  In  the  implicit  framework,  these  have  the  form 

z  =  0,  Ao  e  Ao  C  (3) 


I-AoA 

C 


where  A  £  C"^",  C  G  and  the  structure  Aq  is  still  of  the  form  (1),  but  with  blocks 

which  are  constant,  complex  matrices  rather  than  dynamical  operators.  Aq  could  also  have 
blocks  restricted  to  be  real  as  in  [27],  which  can  be  used  to  capture  real  parametric  uncertainty. 
By  analogy  with  the  dynamic  case,  we  will  say  that  the  implicit  system  (3)  is  stable  if 


Ker 


I- AoA 
C 


=  0  VAo  G  Bao 


(4) 


Condition  (4)  strongly  resembles  the  PBH  test  in  standard  system  theory.  In  fact,  for  the  spe¬ 
cial  case  Ao  =  SI,  stability  is  equivalent  to  detectability  of  the  pair  (C,A).  Further  connections 
will  be  mentioned  in  section  3. 

Different  problems  which  lead  to  a  constant  matrix  formulation  are  now  reviewed. 


2.1  LTI  uncertainty 

Analogously  to  the  standard  input-output  case,  if  the  perturbations  A  are  linear  time- 
invariant  (LTI),  then  robustness  analysis  can  be  reduced  to  a  constant  matrix  test  across 
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frequency  (over  the  unit  circle  or  the  unit  disk).  For  a  finite  dimensional  LTI  map  M, 

will  denote  the  frequency  response,  and  for  causal  M,  M(^),  ^  G  C  will  denote  the  Z-transform 

M(o  =  E*M(t)r. 


Proposition  1  Consider  the  implicit  system  (2),  where  A,  C,  A  are  finite  dimensional  LTI 
maps  in  C{1‘2)-  Let  the  constant  matrix  Aq  have  the  same  spatial  structure  as  A.  Then  system 
(2)  is  robustly  l^-stable  if  and  only  if 


Ker 


-  AoA(e^“) 
C{ei'^) 


0  VAo  6  Bao,  Vw  G  [-7r,7r], 


(5) 


Proposition  2  Consider  the  implicit  system  (2),  where  A,  C,  A  are  causal  finite  dimensional 
LTI  maps  in  jC(l2e)-  Let  the  constant  matrix  Aq  have  the  same  spatial  structure  as  A.  Then 
system  (2)  is  robustly  stable  if  and  only  if 

=  0  VAoGBao,  V|e|<l  (6) 


Ker 


I  -  KoAifi) 

m 


2.2  State  space  representations 


Another  standard  method  to  obtain  a  constant  matrix  formulation  is  by  writing  a  state-space 
realization  for  the  discrete- time  (causal)  maps  A,  (7;  consider  the  joint  realization 


A 

C 


A-x 

Axz 

= 

Azx 

Az 

.  Cx 

<7.  . 

=  X{A^x 

+  Axz 

(7) 


the  state.  By  adding  these  equations  in  implicit  form  to  (2)  we  obtain 


A/  0 
0  A 


4  4 

■^X  ^XZ 

A,^  A, 


(8) 
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which  will  be  represented  by 

with  A5,  Cs  constant  matrices,  and  the  augmented  delay-uncertainty  operator  As  =  diag[X,  A]. 
For  LTI  uncertainty  A,  (9)  can  in  turn  be  reduced  to  a  test  in  terms  of  constant  complex 
perturbations  Asoi  constant  matrix  conditions  in  state  space  can  also  be  given  for  the  LTV 
case.  These  issues  will  be  discussed  in  Section  4. 

2.3  Model  Validation  as  Implicit  Analysis 

Substantial  attention  in  recent  years  has  focused  in  establishing  closer  connections  between 
robust  control  and  system  identification.  In  this  category  fall  the  results  of  [22,  13],  where  a 
model  validation  problem  is  posed  as  a  generalization  of  the  structured  singular  value  /i.  In 
particular,  [22]  considers  a  fi  problem  constrained  to  a  subspace. 

This  is  precisely  the  type  of  extension  provided  by  an  implicit  representation  as  (3), 
as  noted  in  [7].  A  simple  example  is  presented  here  to  illustrate  this  point.  Consider  the 
validation  of  a  linear  regression  model  of  the  form 

y  =  M0  +  d  (10) 

with  ||0||  <  it",  ||d||  <  7.  In  this  equation,  y  is  a.  given  vector  in  R",  and  M  a  given  matrix 
in  R"^™,  both  related  to  experimental  data.  The  validation  problem  is  to  determine  whether 
there  exist  vectors  6  and  d  in  the  allowed  class,  satisfying  (10). 

The  basic  observation  is  that  the  size  constraints  on  6  and  d  can  be  captured  by  implicit 
uncertain  equations  6  =  A.gK  and  d  =  Ad'y,  where  A^,  A^  are  respectively  m  x  1,  n  x  1 
matrices  of  norm  bounded  by  1.  Figure  1  jointly  represents  equations  (10)  and  the  constraints 
by  means  of  an  auxiliary  “input”  w  =  1.  The  existence  of  nontrivial  solutions  for  these 
equations  is  equivalent  to  the  validation  of  the  model  (the  constraint  w  —  1  can  be  obtained 
by  normalization).  This  is  a  stability  question  in  an  implicit  LFT  system,  and  can  be  readily 
reduced  to  the  form  (3)  by  the  method  described  in  Part  1.  For  a  more  complete  discussion 
of  these  issues  in  a  more  general  setting  we  refer  to  [7]. 
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Figure  1:  Validation  of  a  linear  regression  model 


3  A  Structured  Singular  Value  for  Implicit  Systems 

Standard  robust  stability  analysis  for  constant  matrices  is  provided  by  the  structured  singular 
value  [5,  14].  In  the  constrained  case  of  (3),  the  natural  extension  is  given  by  the  following: 

Definition  1  The  structured  singular  value  of  the  implicit  system  (3)  is  defined 

as  follows: 

If  Ker 


min  (Ao)  :  Aq  G  Aq,  Ker 


I-AoA 

C 


=  0  VAo  e  Ao,  define  pii^^{C,A)  =  0,  otherwise 


A  restatement  of  this  definition  is  to  say  that  (3)  is  stable  if  and  only  if  /Uao(C',  A)  <  1. 

Equivalently,  Definition  1  translates  the  analysis  problem  to  the  computation  of  the  func¬ 
tion  iJ,/^g(C,A);  as  in  the  standard  case,  exact  computation  is  in  general  hard  and  one  must 
rely  on  upper  and  lower  bounds.  We  will  only  comment  briefly  here  on  the  lower  bound 
problem,  and  develop  in  detail  the  upper  bound  theory. 

The  lower  bounds  for  the  standard  unconstrained  case  (no  C  equations)  are  based  on  the 
fact  that  fJ,A.iA)  =  maxAgB^(p(AA)),  where  p(-)  denotes  spectral  radius.  Algorithms  which 
resemble  the  power  iteration  for  spectral  radius  have  been  developed  [14,  27],  which  have 
good  performance  on  typical  problems. 
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For  the  constrained  case,  only  eigenvalues  with  eigenvectors  in  the  kernel  of  C  are  relevant. 
In  the  following,  it  will  be  convenient  to  parameterize  this  kernel  by  a  matrix  C]_,  whose 
columns  form  a  basis  for  the  kernel  of  C.  This  leads  to 


Ker 


I- AA 
C 


Ker[Cl  -  AACl]  =  0 


(12) 


Denoting  p^(M,N)  =  max{|/3|  :  /3M  —  Nis  singular}  (  maximum  modulus  of  a  generalized 
eigenvalue  of  M,  iV),  we  have 


p^{C,A)  =  max  p^{C±,AAC_l)  (13) 

These  observations  will  presumably  lead  to  an  extension  of  the  standard  p  lower  bound;  some 
difficulties  arise,  however:  generalized  eigenvalues  do  not  always  exist,  and  also  the  maximum 
in  (13)  need  not  occur  on  the  boundary.  Some  initial  work  is  documented  in  [7]. 

We  will  now  consider  the  upper  bounds  for  this  version  of  the  structured  singular  value. 
The  following  theory  strongly  parallels  that  of  [14]  for  the  standard  case.  We  will  use  the 
same  notation  Y  as  in  Part  I  for  the  real  vector  space  of  hermitian  matrices  which  commute 
with  the  elements  in  Aq,  which  have  the  form  Y  =  diag  [Yi , . . . ,  1^,  UL+iImi )  •  ■  • ,  yL+Flmp]- 
We  recall  that  {Y,Y)  =  J2f=i  tr(yifi)  +  yL+jVL+j  defines  an  inner  product  in  Y.  Also, 
X  and  X  are,  respectively,  the  convex  subsets  of  Y  of  positive  and  nonnegative  scalings. 


Lemma  3  For  fixed  fi  >  0,  the  following  are  equivalent: 

(i)  3X  G  X  :  A*XA  -  (3^X  -  C*C  <  0 
(m)  3X  G  X  :  C^{A*XA-  (3^X)Cl  <  0 


(14) 

(15) 


The  previous  conditions  are  both  Linear  Matrix  Inequalities  (LMIs,  [2])  (strictly  speaking,  (i) 
is  affine  rather  than  linear).  Testing  whether  an  LMI  is  satisfied  is  a  convex  feasibility  problem, 
for  which  interior  point  methods  are  available  [2,  9].  While  version  (i)  is  more  directly  related 
to  robustness  analysis  tests,  (ii)  is  of  lower  dimensionality  and  therefore  preferable  from  a 


7 


computational  point  of  view.  We  define  the  upper  bound  for  jj,, 


AaoCC",  A)  =  inf{;5  >  0  :  (14)  is  satisfied}  (16) 

The  fact  that  Hao(C,  A)  <  jl^^{C,A)  is  a  consequence  of  Theorem  4  below.  We  first  mention 
the  following  remarks: 

•  Combining  the  upper  bound  with  Propositions  1  and  2  provides  tractable  sufficient 
conditions  for  robust  {I2)  stability  in  the  case  of  LTI  perturbations;  for  example, 

3X{uj)  G  X  :  A(e^'“)*X(w)A(e'")  -  X(u;)  -  C'(e''“)*C(e^'“)  <  0  Vw  G  [-x,  x]  (17) 
an  LMI  across  frequency,  guarantees  robust  /2-stability  of  (2)  under  LTI  perturbations. 

•  LMI  (15)  for  =  1  has  appeared  in  previous  work  [10]  on  stabilization  of  input-output 
LFT  systems,  where  it  characterizes  the  so-called  Q-detectability  of  the  pair  {Afi) 
(this  reinforces  the  connection  with  the  PBH  test  mentioned  earlier).  It  is  shown  in 
[10]  that  it  is  equivalent  to  the  existence  of  an  output  injection  matrix  L  such  that 
mixa{X{A  +  LC)X-^)  <  1. 

•  If  the  structure  includes  real  81  blocks  (corresponding  to  parametric  uncertainty),  the 
upper  bound  can  be  improved  in  the  same  manner  as  the  standard  case  (see  [27]). 

For  the  analysis  of  the  upper  bound  we  introduce  a  static  version  of  the  V  set  defined  in  Part 
1.  For  C  G  C”,  let 

$°(C)  =  (AC)i(AC);-cc;  i  =  i...L 
^2+;(C)  =  (AC)U,(AC)r+i  -  J  - 1 . . 

A°(C)  =  diag  [$“(0, . . . ,  $i(0,  ,  •  •  • ,  (18) 

Define  V°  =  {A°(()  :  CC  =  0,  |C|  =  1}  C  Y;  co(V°)  denotes  its  convex  hull. 
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Theorem  4  1.  fj, A)  <12/  and  only  2/ V°  n  X  =  ^. 

2.  The  following  are  equivalent: 

{i)  /2Ao(C',-4)  <  1 

(u)  3X€X:A*XA-X-C*C'<0  (19) 

(222)  co(V°)  n  X  = 

In  this  constant  matrix  case,  the  upper  bound  will  be  strict  in  general;  equivalently,  LMI 
(19)  is  not  a  necessary  test  for  stability.  Referring  to  Part  I,  convexity  of  the  V  sets  played 
an  essential  role  in  the  necessity  results,  and  this  does  not  hold  in  the  static  context:  V°  is 
not  in  general  convex  (so  V°  PI  X  =  ^  does  not  imply  co(V°)  PI  X  =  f). 

In  a  similar  manner  as  in  the  standard  case  [14],  we  now  pose  the  question  as  to  which 
special  A  structures  give  equality  of  jj,  and  f. 

Definition  2  The  structure  A  is  fi-simple  in  the  implicit  case  if  p^(C,  A)  —  p,^{C,A)  for 
any  matrices  A,  C . 

Theorem  5  The  following  structures  are  p-simple  in  the  implicit  case. 

(i)  A  =  :  (i  G  C} 

(ii)  A  =  {diag[lS.i, . . .  A^p]  :  A,  G  with  F  <2,  for  A,C  real. 

(Hi)  A  =  {diag[Ai, . . . ,  Ajr]  :  Aj  G  with  F  <  3. 

In  reference  to  structures  with  only  full  blocks,  the  situation  is  analogous  to  the  standard  case 
of  [14]:  the  bound  is  exact  for  a  maximum  of  3  complex  full  blocks  or  2  real  full  blocks.  The 
only  notable  difference  in  the  implicit  case  is  the  fact  that  the  structure  A  =  {diag[6il,  A2]} 
is  no  longer  /2-simple,  as  shown  in  the  following  example. 
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Example  1  Let  A  =  diag  [^1/2,  Aj],  G  C,  A2  G 


C!  = 


The  top  half  of  Cf  -  AACf  is 


1  0 
0  1 
4  0 
0  4 


1-2^1  0 
0  1-3^1 


1 

'2  O' 

II 

0  3 

3  2 

J 

.2  2_ 

,  so  the  kernel  is  nontrivial  only  for 


b  = 

-  1/2  or  Si 

4 

—  A2 

3 

0 

2 

.  This  can  be  achieved  with  <t(A2)  of  at  least  Ajy/lA. 

A  similar  argument  with  Si  =  1/3,  shows  that  for  a  nontrivial  kernel,  d(A2)  >  V^-  The 
first  perturbation  is  smaller  so  pL^{C,A)  =  \/l3/4  <  1. 


For  the  LMI,  write  X  =  diag[Xo,/2],  with  Xq 


X  y 

y*  z 


.  Some  algebra  gives 


d{x-l)  5(t/  +  2) 
5(r +  2)  8(z-  1) 


C^{A*XA-  X)Cl  = 

For  (20)  to  be  negative  definite,  and  X  >  0,  we  must  have 


24, 


0<a;<l;  0<2<1;  \y\‘^  <  xz-,  \y  +  2\^  < —{1  -  x){l  -  z) 

Zo 


(20) 


(21) 


This  implies  \y\  <  1,  |?/  +  2|  <  1  which  is  impossible,  so  there  is  no  solution  to  LMI  (15)  with 
/3  =  1.  Consequently,  p^{C,A)  <  /1a (C*?^)- 


To  conclude  this  section,  we  relate  the  LMI  test  (19)  to  the  results  of  Part  1.  If  A  and  C 
remain  constant  but  allowed  to  operate  on  I2  signals,  and  Aq  is  substituted  by  an  arbitrary 
structured  operator  on  I2,  Theorem  1,  Part  I  implies  that  (19)  is  a  necessary  test  for  robust 
stability.  This  amounts  to  an  infinite  horizon  augmentation  of  the  constant  matrix  problem 
(3);  for  this  case  where  A  and  C  are  constant,  a  finite  horizon  augmentation  suffices: 

Theorem  6  Let  A,C,  be  constant  matrices,  operating  on  finite  horizon  signals  z  =  {zk}i~i  G 
(C")'^.  Let  A  be  the  class  of  structured  (as  in  {!))  operators  in  (C")'^.  If  d  >  dim(Y),  then 

=  0  VA  G  Ba  ^  (19)  holds  (22) 


Ker 


I-AA 

C 
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The  previous  finite  horizon  augmentation  can  also  be  rewritten  in  matrix  form,  by  defining 
matrices  which  are  d  times  larger  than  A,  C,  and  Aq,  obtained  by  adequate  repetition.  A 
similar  result  (for  standard  //-analysis)  has  been  obtained  using  very  different  methods  in  [1], 
where  the  size  of  the  augmentation  is  n  rather  than  dim{Y)  (these  two  are  the  same  when 
A  consists  of  only  scalar  blocks). 

In  comparison,  in  the  case  where  A  and  C  have  unbounded  memory,  an  infinite  horizon 
augmentation  (to  LTV  operators  on  I2  as  in  Part  I)  is  required. 

4  Analysis  of  state  space  systems 

In  this  section  we  consider  the  state  space  representation  of  (9),  where  AsiCs  are  constant 
matrices,  and  the  structure  As  =  diag[A/  A]  has  a  “special”  first  block,  given  by  the  delay  op¬ 
erator.  As  in  Proposition  2,  we  will  be  dealing  with  causal  operators  A  and  the  corresponding 
notion  of  robust  stability. 

4.1  LTI  Uncertainty 

If  the  uncertainty  A  is  LTI,  robustness  analysis  reduces  once  more  to  a  constant  matrix 
problem,  where  the  corresponding  constant  matrix  structure  As^  has  a  61  first  block. 

Proposition  7  System  (9),  with  A  a  structured  LTI  operator  in  Ccihe);  is  robustly  stable  if 
and  only  if  As)  <  1. 

The  upper  bound  for  //  will  provide  a  computationally  tractable  sufficient  condition  for 
robust  stability  of  the  form 

A*  XsAs  -  X5  -  C*sCs  <  0  (23) 

where  Xs  =  diag[XQ,X]  is  defined  to  commute  with  A5:  Xq  is  a  positive  square  matrix  of 
dimension  equal  to  the  number  of  states,  and  X  G  X.  This  condition  is  in  general  conservative, 
even  in  the  case  where  A  is  unstructured  (one  full  block);  this  is  a  consequence  of  the  fact  that 
in  the  implicit  case,  the  structure  Asg  =  diag[6I,A]  is  not  //-simple,  as  shown  by  Example  1. 
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4.2  LTV  Uncertainty 

Since  (23)  is  a  constant  scales  test  it  can  be  related  to  analysis  with  LTV  perturbations,  as 
considered  in  Part  1.  We  have  the  following: 


Proposition  8  If  (23)  holds,  (9)  is  stable  for  all  As  =  diag[X,A],  A  G  Ba  where  A  is  the 
set  of  arbitrary  causal  bounded  LTV  operators  (A  C  Ccihe))- 


It  is  not  clear  in  general  whether  the  converse  of  Proposition  8  holds,  as  it  does  in  the 
unconstrained  case  (no  0$)',  the  results  of  Part  I  do  not  apply  directly  since  here  the  first 
block  of  As  is  constrained  to  be  the  delay. 

A  special  case  where  the  converse  holds  is  when  the  constraints  Cs  do  not  involve  the 
state  variables  x:  consider  the  state  space  implicit  system  of  Figure  2  a),  which  corresponds 
to  the  case  =  0  in  (8).  The  corresponding  system  (2),  depicted  in  Figure  2  b),  has  static 
C  constraints.  As  was  shown  in  Part  I,  a  general  robust  performance  problem  with  a  finite 
number  of  IQCs  in  the  disturbance  variables  can  be  cast  in  this  special  form. 

Theorem  9  In  reference  to  Figure  2,  let  A  vary  in  the  class  of  structured,  otherwise  arbitrary 
causal  operators  in  Ccihe)-  The  following  are  equivalent: 


(0 

(ii) 

(in) 

(iv) 


(23)  holds  with  As 


A...  A..., 

1 

^zx  -^z 

II 

0  c' 

The  implicit  system  (9)  of  Figure  2  (a)  is  robustly  stable. 

p(As;)  <  1,  and  the  implicit  system  (2)  of  Figure  2  (b)  is  robustly  stable. 

p(A^)  <  1,  and  3A  e  X  :  A(e^“)*XA(e^“)  -  A  -  C*C  <  0  Vcu  G  [-7r,7r] 


(24) 


For  the  case  of  LTV  perturbations  considered  above,  the  same  X  (constant  scaling)  must 
hold  across  frequency,  which  is  a  stronger  condition  than  the  one  given  in  (17)  for  the  LTI 
case.  Solving  (24)  in  a  set  of  frequency  points  gives  a  coupled  LMI  problem.  The  “state-space 
LMI”  (23)  gives  X  in  just  one  LMI,  but  must  handle  an  extra  full  block  in  the  scaling  Xs- 
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(a)  (b) 

Figure  2:  Two  formulations  of  the  robust  stability  problem 

Remark:  Condition  (iv)  can  be  rewritten  (see  the  proof  of  the  theorem)  in  the  form 

X"^ACl{C^XCl)-^  <  1  (25) 

oo 

which  is  especially  adequate  for  synthesis  methods  extending  the  so-called  D  —  K  iteration 
for  ^-synthesis.  In  this  case  A  is  a  function  of  a  controller  K,  and  an  Tioo  synthesis  step  is 
alternated  with  an  analysis  fit  of  X.  These  issues  are  addressed  in  [4]. 

5  Application  to  Robust  7^2  Analysis 

In  this  section  we  apply  the  implicit  analysis  framework  developed  in  this  paper  to  the  problem 
of  analyzing  white  noise  rejection  properties  of  an  uncertain  system. 

For  an  LTI  input-output  system  in  the  absence  of  uncertainty,  the  relevant  measure  is  the 
■^2  norm  of  the  system,  which  arises,  for  example,  as  the  expected  output  power  when  the 
input  is  a  stationary  random  process  with  flat  power  spectrum. 

The  stochastic  paradigm  is  less  attractive,  however,  when  analyzing  systems  subject  to 
additional  sources  of  uncertainty  (parameters,  unmodeled  dynamics)  which  are  usually  ex¬ 
pressed  more  naturally  in  a  deterministic  setting.  Research  addressing  this  “Robust  'H2' 
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problem  [15,  23]  has  faced  the  difficulty  of  analyzing  the  average  effect  of  the  disturbance 
together  with  the  worst-case  effect  of  the  nncertainty.  This  is  the  main  reason  for  which  other 
deterministically  motivated  performance  measures  (Tfoo,  Ti)  have  gained  popularity  in  the 
robust  control  literature. 

For  many  disturbances  arising  in  applications,  however,  a  white  noise  model  is  more 
appropriate  than,  for  instance,  the  class  of  arbitrary  bounded  power  signals  considered  in  the 
Tfoo  problem,  where  the  worst-case  signals  (sinusoids)  are  often  very  unrealistic. 

The  following  treatment  is  based  on  imposing  deterministic  constraints  on  the  disturbance 
inputs  by  considering  set  characterizations  of  white  noise  [16,  17].  These  constraints  can  be 
included  in  an  implicit  analysis  problem  of  the  form  (2),  and  analyzed  within  this  framework. 

5.1  Set  descriptions  for  white  noise 

This  section  shows  how  deterministic  descriptions  of  white  noise  can  be  fit  into  the  implicit 
analysis  framework.  We  will  first  consider  the  case  of  scalar  signals.  In  [17],  sets  of  signals 
with  a  parameterized  degree  of  “whiteness”  are  defined  as 

=  {w  e  I2,  |r„(r)|  <  7r„(0)  r  =  1 . .  .T}  (26) 

where  r^^r)  =  {w,X^w)  is  the  autocorrelation  of  a  scalar  signal  w{t).  Restricting  an  input 
signal  to  such  a  class  is  a  deterministic  method  to  rule  out  highly  correlated  signals  (e.g., 
sinusoids)  which  are  deemed  unrealistic  in  a  particular  problem,  and  which  dramatically 
affect  the  system  gain.  In  fact,  the  worst-case  gain  of  an  input-output  LTI  system  H  under 
signals  in  kFy.T,  denoted  satisfies 

\m\l  <  iiff <  m\l + 7  E  k»(r)i  +  Y, 

t=-T  \'t\>T 

where  r;,(r)  are  the  autocorrelations  of  the  system  impulse  response.  So  for  small  7  (one  could 
use  7  =  0)  and  large  T,  the  worst  case  induced  norm  of  the  system  under  the  autocorrelation 
constraints  approximates  the  7^2  norm.  Alternative  descriptions  in  the  frequency  domain  are 
also  considered  in  [17],  by  restricting  signals  to  have  equal  energy  in  a  number  of  frequency 
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bands,  thereby  enforcing  an  approximately  flat  spectrum.  The  relationship  between  these 
descriptions  and  the  alternative  stochastic  paradigm  is  analyzed  in  [16,  17]. 

The  advantage  of  a  description  such  as  (26)  based  on  a  finite  number  of  quadratic  con¬ 
straints  (as  opposed  to  the  ideal  specification  that  ru,(r)  be  the  delta  function)  was  already 
pointed  out  in  [11].  In  our  setting,  it  corresponds  to  the  fact  that  the  resulting  sets  W-^^t 
can  be  represented  by  the  behavior  of  an  implicit  uncertain  system,  as  was  shown  for  general 
IQCs  in  Part  I. 

We  include  the  following  derivation  to  motivate  the  multivariable  case.  Consider  the 
constraints  in  (26),  for  real  signals: 

±  ^^{t)  <  7  r,„(0)  r  =  1 . .  .r  (28) 

Simple  manipulations  reduce  (28)  to 

2(1  -  7)l|w||^  <  jju;  ±  X^wf  t=1...T  (29) 

For  a  fixed  r  and,  for  instance,  the  plus  sign  in  (29),  let  =  a/2(1  —  7),  =  1  +  A’’.  The 

corresponding  constraint  ||P^u)|[^  <  ||(5+ry||^  is  equivalent  by  Lemma  1,  Part  I  to 

(p;  -  ilQX)w  =  0  (30) 

for  some  contractive  operator  The  same  procedure  can  be  repeated  for  the  minus  sign 
in  (29),  and  for  r  =  1 . .  .T.  The  constraints  (30)  can  then  be  jointly  represented  by  (P  — 
^cQ)w  =  0,  where  P  =  [Pi,P^  .  .P!"]',  Q  =  [Q\,Q^_  . .  .Ql]' ,  and  Ac  =  diag[6X,  6l  . . .  61]. 
The  set  of  white  signals  Wy,T  has  been  represented  in  the  form 

W,,T=  U  KeriP-AcQ)  (31) 

Ac6Ba 

The  previous  construction  can  be  extended  to  the  multivariable  case,  by  considering  the 
autocorrelation  matrix  of  a  vector  valued  signal  w  G  l^,  Rw(t)  —  J2tL-oo  'wit  +  T)w{t)* .  For  w 
to  be  white,  Pw(r)  must  be  0  for  r  7^  0,  and  ^^(0)  must  be  a  multiple  of  the  identity  matrix. 
These  matrix  conditions  could  be  reduced,  entry  by  entry,  to  a  number  of  scalar  constraints, 
and  treated  as  before. 
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For  r  ^  0,  and  using  the  complex  field,  a  simpler  method  is  given  by  operator-valued 
61  blocks  in  the  uncertainty  Let  r  7^  0  be  fixed,  and  define  =  y/2I,  =  (1  i 

Q'^  =  (1  ±  jA’’)/.  Consider  the  following  four  implicit  equations,  where  each  6  is  an  arbitrary 
operator  on  /2- 

(P^  -  611  Ql)w  =  0  (32) 

{P^-6lIQl)w  =  0  (33) 

By  use  of  Lemma  2,  Part  I,  (32-33)  are  equivalent  to 

2\\T]*wf  <  I|7/*(w±  A^w)|p  Vr/GC™  (34) 

2||77*u;|P  <  ||77*(tt?  ±  yA’'ry)|p  E  C""  (35) 

which  in  turn  reduce  to  {r]*  w ,  rj* w)  =  0  Vr/  G  C’”-t=^P„(r)  =  0. 

Comparing  with  the  IQC  formulation,  the  constraints  (32-33)  for  multivariable  white  noise 
correspond  to  matrix-valued  IQCs,  as  mentioned  in  Part  1. 

5.2  Robust  7^2  Performance  Analysis 

The  framework  will  now  be  applied  to  a  problem  of  white  noise  rejection  analysis  in  the 
presence  of  uncertainty.  For  this  purpose,  we  return  to  the  general  setup  of  Section  3.4,  Part 
I,  which  for  convenience  is  represented  in  Figure  3.  We  are  given  an  uncertain  input-output 
system  given  as  an  LFT  between  an  LTI  map  H(X)  and  a  structured  uncertainty  operator  A^. 
The  question  is  to  test  whether  the  worst-case  I2  gain  from  d  to  y  in  the  presence  of  uncertainty 
A„  is  less  than  /?,  when  the  input  signal  d  is  forced  to  satisfy  “whiteness”  constraints  of  the 
form  discussed  in  section  5.1.  These  constraints  are  represented  by  a  map  P  —  AcQ,  on  the 
left  in  the  picture,  where  without  loss  of  generality  P  can  be  chosen  to  be  static.  The  same 
formulation  allows  for  arbitrary  IQCs  applied  to  d.  The  “performance  IQC”  ||r/||  >  /?  ||d||  is 
represented  by  the  block  Ap. 

^  This  provides  a  motivation  for  this  type  of  blocks 
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Figure  3:  Robust  0,2  Analysis 


The  reduced  representation  (2)  for  this  problem  was  already  obtained  in  Part  I.  Adding 
the  scaling  (3  gives 


’  A„ 

0 

0 

'Hn 

^H,2 

O' 

A  = 

0 

Ap 

0 

,  A  = 

H21 

-pH22 

0 

,  C'  = 

0  iF  -/], 

2;  = 

Zp 

0 

0 

Ac 

0 

jQ 

0 

-  . 

(36) 


In  this  case,  the  blocks  Ap  and  Ac  are  already  in  the  class  of  structured  LTV  operators.  We 
will  analyze  the  case  where  A„  is  also  structured  LTV,  which  may  give  a  conservative  answer 


if  it  includes  parametric  or  LTI  uncertainty.  Since  C  is  static,  the  robust  stability  test  is  given 
by  (24),  or  equivalently  (from  Theorem  9)  by  the  state-space  version  (23),  which  is  obtained 
by  writing  state-space  realizations 


'  Ah 

Bhi 

Bh2 

/7(A)  = 

Chi 
_  Ch2 

Dhii 

Dh21 

Dh12 
Dh22  . 

Q(A)  = 


Aq 

Bq] 

Cq 

Dq\ 

The  corresponding  equations  (9)  have  A5  =  diag[\I,  A/,  A^,  Ap,  Ac], 


Ah 

0 

Bhi 

jBh2 

O' 

0 

Aq 

0 

PQ 

0 

Chi 

0 

Dhii 

^Dh12 

0 

(37) 

Ch2 

0 

Dh12 

jDh22 

0 

0 

0 

jDq 

OJ 

0 

0 

0 

-P  - 

V 

(38) 
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Let  /3opt  be  the  infimum  of  the  values  of  the  parameter  ^  such  that  LMI  (23)  is  feasible;  this 
is  a  measure  of  the  worst-case  gain  under  uncertainty  and  autocorrelation  constraints. 
Asymptotically,  as  the  number  of  constraints  increases,  the  process  converges  down  to  a 
robust  7^2  performance  measure,  so  that  a  finite  number  of  constraints  always  gives  an  upper 
bound. 

5.3  Examples 

We  will  present  two  simple  examples  to  demonstrate  the  machinery,  applied  to  problems 
involving  the  'H2  norm. 

5.3.1  An  example  without  uncertainty 

The  first  example  consists  of  calculating  the  7^2  norm  of  the  transfer  function  if  (A)  = 
using  this  approach.  There  are  of  course  exact  ways  to  compute  the  H2  norm,  which  give  a 
result  of  1/a/3  =  0.577;  this  example  is  included  for  verification  purposes. 

The  process  described  above  was  performed  with  a  number  T  of  autocorrelation  constraints 
(for  7  =  0).  The  feasibility  of  LMI  (23)  was  checked  using  the  software  package  LMI-Lab  [9]. 
Figure  4  depicts  (3 opt  as  a  function  of  T.  Starting  at  T  =  0  with  the  unconstrained  (Tioo)  norm 
which  is  1,  (3 opt  asymptotically  converges  to  the  7^2  norm,  as  expected,  at  a  rate  consistent 
with  bound  (27). 
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5.3.2  Robust  H2  example 


We  consider  the  standard  SISO  feedback  system  of  Figure  5,  where  the  plant  P  is  subject  to 
multiplicative  uncertainty.  We  wish  to  analyze  sensitivity  of  the  tracking  error  e  to  a  white 
noise  disturbance  appearing  in  d  (which  could  be  due,  for  example,  to  sensor  noise).  The 
map  from  d  to  e  (sensitivity  function)  of  the  uncertain  system  is  given  by 

9  =  _ i _  ('39') 

1  +  Pir(l  +  WA„)  l  +  WToA^  ^  > 

where  So  =  y+pk  “  ipPK  nominal  sensitivity  and  complementary  sensitivity 

functions. 


Figure  5:  Rejection  of  sensor  noise 


LTI  uncertainty: 

If  the  perturbation  A^  is  assumed  be  LTI,  for  this  simple  case  the  worst-case  7^2  norm 
can  be  computed  in  the  frequency  domain;  this  will  allow  us  to  evaluate  the  results  obtained 
from  the  analysis  framework.  Assuming  that  ||kFro||^  <  1,  we  have 

max  |5(e-’“)|  =  - — }  .  . .  (40) 

This  fact  has  been  typically  used  (see  [6])  to  show  that  the  worst  case  Hoo  norm  of  the  system 


is  given  by 


l-IWTol 


.  Here  we  will  use  it  to  obtain  a  worst  case  Ho  norm  of 


l-\WTo\ 


;  for 


this  we  allow  A^  to  be  a  noncausal  (£00)  operator,  which  can  achieve  bound  (40)  for  every 
frequency. 

We  now  choose  K  =  2,  P  =  and  W  =  0.25.  These  values  were  chosen  so  that 

the  uncertainty  affects  the  sensitivity  in  a  significant  way;  this  is  exhibited  in  Figure  6,  where 
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the  lower  curve  indicates  the  nominal  sensitivity  function,  and  the  upper  curve  the  worst-case 
sensitivity  from  (40).  We  obtain  the  values 


5o 

l-\WTo\  ^ 
So 

1  -  iWTol  , 


6.50 

2.39 


(41) 

(42) 


Figure  6:  Nominal  and  worst-case  sensitivity  functions  (magnitude) 

LTV  uncertainty 

Exact  analysis  for  an  arbitrary  LTV  operator  can  be  obtained  from  the  procedure 
described  in  Section  5.2.  Figure  (7)  shows  the  corresponding  plot  of  f3opt  (obtained  using 
LMI-Lab)  as  a  function  of  the  number  of  correlation  constraints  T. 


20 


For  r  =  0  (no  constraints)  we  retrieve  the  value  from  (41)  for  the  worst  case  Jico  norm 
(it  is  well  known  that  in  this  unstructured  case,  the  worst-case  perturbation  is  LTI).  As  T 
increases,  we  approach  the  worst-case  gain  under  white  noise  signals.  The  asymptotic  value 
also  appears  to  coincide  in  this  case,  with  the  value  (42)  obtained  from  LTI  uncertainty,  which 
is  plotted  for  comparison 

Although  for  this  case  the  frequency  domain  method  is  much  simpler,  it  does  not  generalize 
to  multivariable  systems  or  to  structured  uncertainty.  The  procedure  based  on  the  implicit 
framework  applies  in  principle  to  any  case,  although  for  this  method  to  be  practical  in  large 
problems,  improvements  in  the  efficiency  of  LMI  solvers  are  required. 

6  Conclusions 

Implicit  representations  have  been  shown  to  be  an  attractive  general  framework  where  var¬ 
ious  forms  of  system  uncertainty,  performance  requirements  and  signal  constraints  can  be 
expressed.  A  robustness  analysis  theory  has  been  developed  which  includes  what  is  available 
for  the  standard  input/output  setting,  and  enhances  its  domain  of  applications  to  include 
overconstrained  problems,  exemplified  in  this  paper  by  robust  7^2  performance  analysis. 

The  computational  properties  of  this  extension  are  similar  to  those  of  the  standard  case. 
Conditions  obtained  in  terms  of  LMIs  lead  in  principle  to  tractable  computation,  but  the 
size  of  the  problems  is  also  a  concern.  In  this  respect,  tests  with  LTV  uncertainty  yield 
either  coupled  LMIs  across  frequency  such  as  (24),  or  a  large  full  block  of  multipliers  for  the 
state  space  version  (23).  Additional  research  is  required  on  practical  methods  for  these  large 
problems,  and  also  for  the  case  of  mixed  LTI/LTV  uncertainty,  which  arises  naturally  in  this 
setting  (for  example,  when  in  (36)  is  LTI).  Here  the  coupling  is  not  easily  avoided;  some 
initial  work  on  this  problem  is  reported  in  [20]. 

One  of  the  main  reasons  to  adopt  the  LFT  framework  is  that  it  allows  for  the  consideration 
of  highly  structured  (e.g.  real  parametric)  uncertainty  which  is  not  captured  by  IQCs.  For 

^This  is  not  a  general  fact;  other  examples  (see  [18])  exhibit  a  gap  between  LTI  and  LTV  uncertainty. 
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these  cases  the  exact  analysis  conditions  in  terms  of  //  will  have  increased  computational 
complexity,  as  in  the  standard  case  where  they  are  known  to  be  NP  hard  [3].  Although  this 
implies  unacceptable  computation  time  in  the  worst  case,  p  lower  bounds  [14,  27]  have  proven 
to  be  efficient  on  “typical”  problems.  Their  extension  to  the  implicit  framework  is  a  direction 
of  future  research. 

If  successful,  this  extension  can  have  an  enormous  impact  in  problems  involving  data, 
such  as  the  model  validation  problem  mentioned  in  Section  2.3,  and  the  corresponding  system 
identification  problem  (see  [7]).  The  implicit  LFT  framework  would  then  appear  as  the  natural 
setting  for  unifying  modeling,  analysis,  model  validation  and  system  identification,  under  a 
common  set  of  mathematical  and  computational  tools. 

Another  natural  extension  of  the  results  in  this  paper  is  the  question  of  synthesis  of 
controllers  for  robust  performance  in  this  setting.  In  [4],  it  is  shown  how  standard  “D-K 
iteration”  methods  for  /^-synthesis  extend  to  this  formulation,  and  allow  in  particular  for 
design  of  controllers  for  Robust  0,2  performance. 
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Appendix:  Proofs 

Propositions  1  and  2 

The  proof  follows  easily  from  Proposition  5  in  Part  I.  We  sketch  the  argument  for  the  case  of 
Proposition  1,  the  other  one  is  analogous.  For  any  finite  dimensional  LTI  A,  A(e-^“')  G  Bao 
for  every  a;;  if  condition  (5)  holds, 

-  A(e^'‘")A(e^'“) 

C{ei^) 


(43) 
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is  full  column  rank  for  all  w  G  [— 7r,7r],  therefore  I2  stability  follows  from  Prop.  5,  Part  I. 
Conversely  if  (5)  fails  at  some  uq,  Aq,  it  is  easy  to  construct  an  LTI  perturbation  A  G  Ba 
such  that  A(e'^‘*''’)  =  Aq,  which  violates  stability  from  Prop.  5,  Part  I. 

Lemma  3 

Since  CC\_  =  0,  (i)  implies  (ii).  If  (ii)  holds,  there  exists  X  G  X  such  that  A*XA  —  X  <  0 
on  the  kernel  of  C.  By  continuity,  there  exists  e  >  0  such  that  {{A*XA  —  X)u,  u)  <  0  for  all 

Ibll  =  lic^ll^  <  e- 

Now  choose  77  >  0  such  that  Xmax{A.*T]XA  —  rjX)  <  e. 

This  gives  {{A*r]XA  —  rjX  —  C*C)v,  n)  <  0  for  all  n  7^  0,  so  rjX  solves  (i). 

Theorem  4 

1.  Assume  iJ,^^{C,A)  >  1.  Then  there  exists  Aq  G  Baq  such  that  (3)  has  a  nontrivial 
kernel;  let  (  of  norm  1  be  in  the  kernel.  Then  CC  =  0,  AqA^  =  C;  focusing  on  the  blocks 
of  Aq,  we  obtain 

6J(A0i  =  Ci  =>  {AOiiAO*>QQ  $“(O>0  i  =  l...L 
Al+j{A0l+j  =  a+j  ||(^C)n+ill  >  lia+ill  ^  <+j(0  >0  j  =  l...F  ^  ^ 

Therefore  the  matrix  A'’(C)  is  in  D  X.  The  converse  follows  similarly. 

2.  The  equivalence  of  (i)  and  (ii)  is  obvious  from  the  definition  of  fiAo(C,  A). 

Let  X  >  0  solve  (ii).  For  any  (  G  C”,  =  0,  some  algebra  shows  that 

{X,A\0)  =  Etr(Vi<S?(C))  =  C(A-XA  ^X-  C'CX  <  0  (45) 

i=l  j=l 

Also,  (X,  y)  >  0  for  all  y  G  X.  Therefore  the  hyperplane  (X,  y)  =  0  in  Y  separates  the 
sets  V°  and  X,  which  implies  their  respective  convex  hulls  co(V°)  and  X  are  disjoint, 
proving  (Hi).  Conversely,  if  co(V‘’),  X  are  disjoint,  a  separating  hyperplane  can  be 
found  leading  back  to  (ii). 
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A  Lemma  from  Convex  Analysis 

Lemma  10  Let  IC  C  V ,  where  V  is  a  d  dimensional  real  vector  space.  Every  point  in  co(K) 
is  a  convex  combination  of  at  most  d  +  1  points  in  K;  for  K  compact,  every  point  in  the 
boundary  of  co{K)  is  a  convex  combination  of  at  most  d  points  in  K . 

The  first  statement  is  a  classical  result  from  Caratheodory  (see  [21]),  which  implies  that  for 
every  v  G  co(K),  there  exists  a  simplex  of  the  form 

{rf+l  d+l 

X]  :  a*  >  0,  X)  ^ 

k=l  k=l 

with  vertices  v^  G  K,  which  contains  v.  If  the  Vk  are  in  a  lower  dimensional  hyperplane,  then 
d  points  will  suffice.  If  not,  then  every  point  in  S{vi, . . . ,  Uai+i)  corresponding  to  Oj,  >  0 
will  be  interior  to  S{vi, . .  .,nd+i)  C  co{K).  Therefore  for  points  v  in  the  boundary  of  co{K), 
one  of  the  akS  must  be  0  and  a  convex  combination  of  d  points  will  suffice. 

Theorem  5 

(i)  In  the  case  A  =  61,  if  pi£i.{C,  A)  <  1  then  (C,  ^4)  is  detectable  in  the  usual  system  theoretic 
sense,  so  there  exists  an  output  injection  L  such  that  p{A  +  LC)  <  1.  From  Lyapunov  theory 
this  implies  there  exists  X  >  0  such  that 

{AA  LCyX{AE  LC)- X  <Q  (46) 

Multiplying  on  the  left  and  right  by  Cx,  Cf  gives  C^{A*XA  —  X)Cf  <  0  which  implies 
Aa(C,  A)  <  1. 

(ii)  The  only  nontrivial  case  is  T  =  2.  Let  A,  C,  X  =  diag[Xi,  X2]  be  real  matrices.  To 
analyze  this  case  we  must  consider  a  real  version  of  the  V°  set,  of  the  same  form  as  (18)  but 
with  C  G  R”.  Consider  the  n  x  r  matrix  F  =  CX  parameterizing  the  kernel  of  C  (C  =  Tu), 
and  assume  it  is  isometric.  Then  V°  can  be  rewritten  as 

V°  =  {A(?;)  =  [CTi(n),  (72{v)],  v  G  R^  ||u||  =  1}  (47) 
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where  crj(v)  =  v'ffjv,  and  Hj  =  (yir)j(Ar)j  — (r)j(r)j  is  areal,  symmetric  matrix  for  j  =  1,2. 
To  prove  that  this  structure  is  ^-simple  is  equivalent,  by  Theorem  4  to  the  fact  that 

n  X  =  0  =>  co(V“)  n  X  =  0  (48) 


for  any  A,  C,  or  equivalently  for  any  symmetric  Hi,  H2.  In  this  case  X  =  the  closed 

first  quadrant  in  R^.  We  have  therefore  restated  the  problem  as  a  geometric  condition  on 
the  range  of  two  real  quadratic  forms.  In  this  notation,  (48)  is  equivalent  to  an  “S-procedure 
losslessness”  theorem  from  Yakubovich  [25];  since  this  literature  is  not  easily  accessed  we 
include  a  proof  which  is  based  on  some  modifications  to  the  parallel  results  of  [5]. 

Let  P  —  A(up),  Q  —  K.{vq),  be  two  distinct  points  in  V°  {vp,Vq  G  R”,  ||up||  =  HugH  =  1). 
Define 


M  = 

v'p 

[vp  Vq 

1 

.  . 

i 

II 

V'p 

Hi 

\  Vp  Vq 

1 

l^Q. 

L  J 

E  =  {[  rj'Hirj  ,  t]'H2T)  ],  ??  G  R^  ||7?||  =  l} 


(49) 

(50) 

(51) 


Then 


•  E  C  V“.  This  follows  from  the  fact  that  if  ||77||  =  1, 


Vp  Vq 


M  ^j]  =1  from  (49). 


P,Q  £  E.  For  P  set  r/p  =  M ^ 
and  analogously  for  Q . 


,  which  verifies  ||7?p||^  = 


1  0 


M 


=  Ihnll  =  1, 


•  FI  is  an  ellipse  in  R^  (which  may  degenerate  to  a  segment) 
Parameterize  t]  =  {cos{d),sin{0)),  6  G  [— 7r,7r].  If  Hj  = 


cij  bj 


bj  Cj 


then 


vHjV=^-^ 


cos{26) 

[2  j 

sin(29) 

J  =  l,2 


(52) 


This  implies  that  E  is  the  image  of  the  unit  circle  by  an  affine  map,  an  ellipse. 
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We  have  shown  that  given  two  points  in  V°,  there  exists  an  ellipse  E  C  through  those 
points.  Now  we  return  to  (48).  If  co(V°)  D  (R+)^  ^  0,  since  V°  is  bounded  and  (R+)^  is  a 
cone,  there  exists  a  point  in  the  boundary  of  co(V“)  which  falls  in  the  first  quadrant.  Using 
Lemma  10  there  exist  two  points  P,Q  in  V°  such  that  the  segment  PQ  intersects  the  first 
quadrant.  But  then  the  corresponding  ellipse  E  will  intersect  (R"^)^;  (the  geometric  picture 
is  given  in  Figure  8(ii)).  This  implies  fl  (R+)^  /  0. 


P 


Figure  8:  Illustration  to  the  proofs 


(iii).  We  consider  the  case  F  =  3,  the  others  follow  similarly.  The  same  procedure  as  in  (ii) 
yields 

V“  =  {A(u)  =  [ai{v),a2(v),a3(v)]  eR^,v  £  C ,  ||u||  =  1}  (53) 

where  (Tj(v)  =  v*HjV,  and  Hj  are  complex  hermitian  forms  in  C’’.  Similarly,  we  must  show 
the  geometric  result 

V°n(R+)^  =  0=^co(V°)n(R+)^  =  0  (54) 

for  any  Hi,  H2,  H3.  Once  again,  this  result  appears  in  the  “S-procedure”  formulation  [8]. 
The  following  proof  is  based  on  [5].  In  particular,  it  is  shown  in  [5]  (analogously  to  (52)) 
that  for  the  case  r  =  2,  the  set  V°  is  the  image  of  the  unit  sphere  in  R^  by  an  affine  map 
g  :  R^-^R®.  This  gives  an  ellipsoid  E  (with  no  interior)  in  R®,  which  could  also  degenerate  to 
a  projection  of  such  an  ellipsoid  in  a  lower  dimensional  subspace. 

Given  two  distinct  points  P  =  A{vp),  Q  =  A{vq)  in  V°,  an  analogous  construction  as 
the  one  given  in  (49-51)  (with  analogous  proof)  shows  that  there  is  such  an  ellipsoid  i?  C  V° 
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through  the  two  points. 

Assume  now  that  co(V°)  fl  (R+)^  ^  0.  Picking  a  point  in  the  the  boundary  of  co(V°), 
Lemma  10  implies  that  there  are  3  points  P,Q,R  in  such  that  some  convex  combination 
S  =  aP  +  l3Q  +  jR  falls  in  (R+)^.  Geometrically,  the  triangle  PQR  intersects  the  positive 
“octant”  at  S. 

Claim:  S  lies  in  a  segment  between  2  points  in  V°. 

This  is  obvious  if  P,Q,  R  are  aligned  or  if  any  of  a,f3,  7  is  0.  If  not,  consider  the  following 
reasoning,  illustrated  in  Figure  8.  Write 

S  =  aP  +  (3Q  +  jR  =  aP  +  {f3  +  l) +  jR)  =  aP  +  (/?  +  'y)T  (55) 

where  T  lies  in  the  segment  QR.  Now  consider  the  ellipsoid  P  C  V°  through  Q  and  R.  If 
it  degenerates  to  1  or  2  dimensions,  then  T  E  E  C  and  the  claim  is  proved.  If  not,  T  is 
interior  to  the  ellipsoid  E.  The  half  line  starting  at  P,  through  T  must  “exit”  the  ellipsoid  a 
point  U  E  E  CV°  such  that  T  is  in  the  segment  PU .  Therefore  S  in  the  segment  PP,  and 
P,  P  e  V°,  proving  the  claim. 

To  finish  the  proof,  we  have  found  two  points  in  V°  such  that  the  segment  between  them 
intersects  (R+)^  .  The  corresponding  ellipsoid  E  C  between  these  points  must  clearly  also 
intersect  (R+)^.  Therefore  V°  fl  (R+)^  7^  0. 

Theorem  6 

Sufficiency  of  condition  (19)  is  a  consequence  of  the  more  general  result  given  in  Part  I;  the 
same  argument  as  exhibited  in  the  proof  of  Theorem  1  applies  to  this  finite  horizon  setup. 

For  the  necessity,  assume  (19)  does  not  hold.  Therefore  co(V°)nX  7^  0,  and  we  can  choose 
a  point  in  the  boundary  of  co(V°),  which  belongs  to  X.  Since  V°  is  in  a  d  dimensional  real 
vector  space  Y,  Lemma  10  implies  that  there  exists  d  points  in  V°  whose  convex  combination 

d 

ockMCk)  =  Xo>0,  (56) 

k  =  l 
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with  (k  e  C",  C(k  =  0,  Q!jt  >  0,  and  J2k  “fc  =  1-  Define  ^  =  {zk}i=i-,  Zk  =  y/^Ck-  Since  A,  C 
are  static,  we  find  that  Cz  =  0,  and 


J2(Azk)iiAzk)*  -  (zkUzk)*  >  0 


fc=i 


\\{Az),^j\f  >  \\iz),^j\f 


i  =  1 . .  .L 


j  =  1...F 


(57) 

(68) 


Now  we  apply  Lemmas  1,  2  from  Part  I  to  conclude  there  exists  A  G  Ba,  structured  operator 

"  I- AA^ 

C 


in  (C”)‘^,  such  that  AAz  =  Therefore 


z  =  0  violating  the  hypothesis. 


Propositions  7  and  8 

Proposition  7  follows  by  the  same  arguments  as  Proposition  2. 

Proposition  8:  Since  Cs  is  static,  Theorem  12,  Part  I  is  in  force  and  LMI  (23)  implies  that 
has  a  left  inverse  in  Ccihe)  for  any  structured  As  G  jCcihe),  ll^sH  <  1?  thus  in 
particular  for  As  =  diag[XI ,  A],  with  A  G  Fcihe)-,  ||A||  <  1. 


I  -  As  A 
C 


Theorem  9 


(i)  =>  (m)  This  is  a  special  case  of  Proposition  8. 

I  -  AsAs 
Cs 

As  =  diag[AI,  A],  A  G  Ba-  Setting  As  =  diag[AI,  0]  implies  that  (/—  A^^,)  has  a  left  inverse 


(ii)  =>  (in)  From  (m). 


has  a  left  inverse  in  Cc(,he)  (causal,  finite  gain)  for  every 


r  n 

I  0  o' 

r  7  -  AA. 

0  1 

I  —  AsAs 
Cs 

= 

-AA,,(/- AA,)-i  /  0 

0  /-AA(A) 

0  0  / 

0 

c 

in  Ccihe),  so  p{Ax)  <  1.  Now  for  any  fixed  As,  the  identity 

I  -{I-AA,)-^AA,, 
0  I 

(59) 

implies  that  the  second  term  in  the  right  hand  side  of  (59)  has  a  left  inverse  in  Cc{he)-  So 
I  —  AA 

has  a  left  inverse  in  Cc{ke)- 

o 

(in)  =>  (iv)  This  is  a  direct  application  of  Theorem  11,  Part  I. 

(iv)  (i)  Let  the  columns  of  form  a  basis  for  the  kernel  of  C;  (iv)  leads  to 


CxA(F‘^)*XA(e^‘^)Cl  <  CxXCl  Vw  G  [-7r,7r] 


(60) 
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Since  C±XC\_  >  0,  pre  and  post  multiplication  by  (Cj.XC'X)”^  gives  [jA(A)||^  <  1,  where 
^(A)  =  AC]_{C±XC*^)~^ ,  as  remarked  in  (25).  We  will  now  change  notation,  redefining 

C]_{C±XC\)~^  as  C*^,  since  it  still  spans  the  same  column  space.  So  Cj_XC\  =  I,  and 
i4(A)  =  X^AC]_,  which  has  state  space  realization 


Ax 

■^xz 

Ax 

Ax,  Cl 

■^ZX 

A 

XU,x 

X^ACl 

(61) 


It  is  well  known  that  p{A^)  <  1,  ||^(A)||^  <  1  implies  the  existence  of  a  solution  Xq  >  0  to 


A  A  ' 

■^x  -^xz 

* 

'Xo  o' 

A  A 

•^x  -^xz 

o 

o 

1 _ 

■^zx  ^z 

o 

1 _ 

A^x 

L  J 

- 1 

O 

(62) 


Substituting  the  expressions  for  ,  and  using  C^XC*^  =  /,  (62)  leads  to 


I  0 
0  Cx 


Xo  0 

0  X 


As  — 


Xo  0  1\ 

0  x\j 


I  0 
0  Cl 


<  0 


Since  Cs  =  [0  C],  Cs±  = 


I  0 
0  Cx 


,  so  setting  Xs  =  diag[Xo  X]  gives 


Csi.{A*sXsAs  -  X5)C^x  <  0 


(63) 


which  implies  (i)  from  Lemma  3. 


References 

[1]  Bercovici  H.,  Foias  C.,  Tannenbaum  A.,  Structured  Interpolation  Theory,  Operator  The¬ 
ory  Adv.  and  App.,  47,  pp;195-220,  1990. 

[2]  Boyd,  S.P.,  L  El  Ghaoui,  E.  Eeron,  and  V.  Balakrishnan,  Linear  Matrix  Inequalities  in 
System  and  Control  Theory,  SIAM,  Philadelphia,  Pennsylvania,  1994. 

[3]  Braatz,  R.,  Young,  P.,  Doyle,  J.,  Morari,  M.,  Computational  complexity  of  p  calculation, 
IEEE  TAG  May,  1993. 

[4]  D’Andrea  R.,  Paganini  F.,  Controller  Synthesis  for  Implicitly  Defined  Uncertain  Systems 
Proc.  1994  CDC,  Orlando,  Florida.,  pp.  3679-3684. 


29 


[5]  Doyle,  J.,  Analysis  of  feedback  systems  with  structured  uncertainty,  lEE  Proceedings, 
129,  242-250. 

[6]  Doyle  J.,  Francis  B.A.,  Tannembaum  A.,  Feedback  Control  Theory,  Macmillan,  1992 

[7]  Doyle  J.,  Newlin  M.,  Paganini  F.,  Tierno  J.,  Unifying  Robustness  Analysis  and  System 
ID,  Proc.  1994  CDC,  Orlando,  Florida,  pp.  3667-3672 

[8]  Fradkov  A.,  Yakubovich  V.A.,  The  S-procedure  and  duality  theorems  for  nonconvex  prob¬ 
lems  of  quadratic  programming,  Vestnik.  Leningrad.  Univ.  (1)  1973,  81-87  (in  Russian). 

[9]  Cabinet,  P.,  LMI  Lab,  A  Package  for  Manipulating  and  Solving  LMI’s,  May  1993. 

[10]  Lu  W.M.,  Zhou  K.,  Doyle  J.C.,  Stabilization  of  LFT  Systems,  Proceedings  1991  CDC, 
Brighton,  England,  pp.  1239-1244. 

[11]  Megretski  A.,  S-Procedure  in  Optimal  Non- Stochastic  Filtering,  Tech.  Rep.  Code 
Trita/Mat-92-0015,  Royal  Inst,  of  Tech.,  Sweden. 

[12]  Megretski  A.,  Treil  S.,  Power  Distribution  Inequalities  in  Optimization  and  Robusntess 
of  Uncertain  Systems,  Journal  of  Mathematical  Systems,  Estimation  and  Control,  Vol  3, 
No.3,  pp  301-319,  1993. 

[13]  Newlin,  M.,  Smith,  R.,  A  Generalization  of  the  Structured  Singular  Value  and  its  Appli¬ 
cation  to  Model  Validation  to  appear  IEEE  Trans.  Automat.  Control. 

[14]  Packard  A.,  Doyle  J.C.,  The  Complex  Structured  Singular  Value,  Automatica,  Vol.  29, 
No.  1,  pp.  71-109,  1993. 

[15]  Packard  A.,  Doyle  J.C.,  Robust  Control  with  an  Performance  Objective,  Proceedings 
1987  ACC,  Minneapolis,  MN,  pp  2141-2145. 

[16]  Paganini  F.,  Set  Descriptions  of  White  Noise  and  Worst  Case  Induced  Norms,  Proceed¬ 
ings  1993  CDC,  San  Antonio,  Texas,  pp.  3658-3663. 


30 


[17]  Paganini  F.,  White  Noise  Rejection  in  a  Deterministic  Setting,  CDS  Tech.  Rep.  CDS94- 
015,  California  Institute  of  Technology,  1994,  submitted  to  IEEE  Trans  A.C. 

[18]  Paganini  F.,  D’Andrea  R.,  and  Doyle  J.,  Behavioral  Approach  to  Robustness  Analysis 
Proceedings  1994  ACC,  Baltimore,  MD.,  pp.  2782-2786. 

[19]  Paganini  F.,  Doyle  J.,  Analysis  of  Implicitly  Defined  Systems  Proceedings  1994  CDC, 
Orlando,  FL.,  pp.  3673-3678. 

[20]  Paganini  F.,  Analysis  of  Systems  with  Combined  Time  Invariant/Time  Varying  Struc¬ 
tured  Uncertainty,  submitted  to  1995  ACC. 

[21]  Rockafellar  R.  T.,  Convex  Analysis,  Princeton  University  Press,  Princeton,  New  Jersey, 
1970. 

[22]  Smith,  R.,  Doyle,  J.,  Model  validation:  a  connection  between  robust  control  and  identifi¬ 
cation  ,  IEEE  TAC,  July,  1992. 

[23]  Stoorvogel  A.A.,  The  Robust  TL2  Control  Problem:  A  Worst-Case  Design,  IEEE  Trans. 
A.C.,  Vol  38,  9,  pp.  1358-1370,  1993. 

[24]  Willems  J.C.,  Paradigms  and  Puzzles  in  the  Theory  of  Dynamical  Systems,  IEEE  Trans. 
A.C.,  Vol.  36,  pp.  259-294,  1991. 

[25]  Yakubovich,  V.A.,  S-procedure  in  nonlinear  control  theory,  Vestnik  Leningrad  University: 
Mathematics,  pp.  62-77,  1971,  (in  Russian) 

[26]  Yakubovich,  V.A.,  Dichotomy  and  absolute  stability  of  nonlinear  systems  with  periodically 
nonstationary  linear  part.  Systems  and  Control  Letters,  11(3),  pp.  221-228,  1988. 

[27]  Young  P.,  Robustness  with  Parametric  and  Dynamic  Uncertainty  PhD  thesis,  California 
Institute  of  Technology,  1993. 


31 


